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We present a simple off-lattice hard-disc model that exhibits glassy dynamics. The inherent 
structures are enumerated exactly, transitions between metabasins are well understood, and the 
particle configurations that act to facilitate dynamics are easily identified. The model readily maps 
to a coarse grained dynamic facilitation description. 
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When a liquid is cooled fast enough to avoid crystal- 
lization, it becomes a supercooled fluid before eventu- 
ally forming an amorphous solid. The feature of this 
transformation that is least understood is the lack of any 
structural signature associated with the dramatic slow- 
ing down of relaxation phenomena that occurs just prior 
to the temperature where the solid structure is frozen 
in. While there are many theoretical models describing 
glassy dynamics, three important concepts include the 
landscape paradigm, dynamic heterogeneity and kinetic 
facilitation. 

The potential energy landscape (PEL), which was orig- 
inally introduced by GoldsteinP] and later formalised 
by Stillinger and Weber |J| in terms of inherent struc- 
tures, describes the total system as a single point mov- 
ing through the high-dimensional iV-body potential en- 
ergy function of the configurational coordinates. While 
the topography of this surface can be essentially charac- 
terized by local minima connected by saddle points, the 
exponential increase in the number of distinct minima 
with system size means that a complete description of 
the PEL is only possible for small systems Neverthe- 
less, this approach has become an important theoretical 
and computational tool for the study of glasses because 
much of the interesting supercooled and glassy behav- 
ior can be connected to the statistical properties of the 
landscape. Some recent examples include the investiga- 
tion of thermodynamic finite size effects 0; the existence 
of an ideal thermodynamic glass transition at positive 
temperature^, |^ and the connection between the con- 
figurational entropy, or number of minima accessible to 
the fluid, and structural relaxation!^, via the Adam and 
Gibbs relation 8] . 

Experiments and simulations have shown that spa- 
tially heterogeneous dynamics is a general phenomenon 
of supercooled fluids nearing the glass transition^. In 
particular, simulations show there are large domains 
of particles that move very slowly while other more 
mobile particles tend to cluster together in string-like 
arrangements [Toj . This raises an important question; 



How are these dynamic features connected to the con- 
figurational structure of the fluid? Studies of the propen- 
sity of a particle in a single configuration to move a long 
distance have confirmed the presence of a definite config- 
urational component to dynamical heterogeneity llj, but 
they also show that this is not connected to any obvious 
structural indicators like free volume. 

A recent theoretical development IT^ that focuses 
on the important role of dynamic heterogeneity in glassy 
systems is based on dynamic facilitation^^. The cen- 
tral idea is that the dynamical structure observed in su- 
percooled liquids is the consequence of local dynamical 
rules that restrict trajectory space, rather than any par- 
ticular property of the static interaction potential. To 
exemplify this, Garrahan and Chandler have used the 
Fredickson-Andersen 14J and EastJ^] spin lattice mod- 
els, which only allow spins the opportunity to flip if cer- 
tain local configurational constraints are statisfied, to ex- 
amine some generic features like "hopping time" distribu- 
tion functions 13] and exchange times 16] that might be 
measured in experiment. However, these coarse grained 
models avoid the question concerning the nature of the 
configurational component that gives rise to the local dy- 
namical rules in structural glasses. 

In this letter we present a simple off-lattice model for 
which we can clearly identify the local configurational 
component and dynamical rules that give rise to dy- 
namic heterogeneity and facilitation. We also obtain a 
complete description of the landscape for the model, in- 
cluding a complete enumeration of all the inherent struc- 
tures, how they are simply connected and organised into 
metabasins, making this a system for which the main 
glassy phenomenologies, the potential energy landscape, 
dynamic heterogeneities and kinetic facilitation can be 
investigated at the same time. We explore how these ap- 
proaches describe the glassy behavior of the model and 
how, in this simple case, all three are interconnected. 

The model consists of N 2-dimensional hard discs of 
diameter cr confined between two hard walls (lines) of 
length L separated by a distance H = [h + 1) where 
lengths are given in units of a. The particle-particle and 
particle-wall interaction potentials are given by 
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respectively, where rij = |rj — ri| is the distance between 
particles and ry is the component of the position vector 
for a particle perpendicular to the wall. The occupied 
volume is then Z2d = Nttct'^ /{4:L{h + 1)), but because of 
the quasi 1-dimensional nature of system, it is useful to 
use the occupied length z = Ncr/L = Z2d^[h + l)/iTa. 
Finally, by restricting the channel diameter to h < ■\/3/4 
we ensure that only nearest neighbors can interact. This 
also prevents the particles from passing each other. We 
perform event driven molecular dynamics simulations 
on a system with h — 0.866. Time t has units of 
CT(m/A;T)i/2. 

To construct the hard particle equivalent to the poten- 
tial energy landscape we need to obtain the distribution 
of jammed packed configurations with respect to their 
packing density 0,^3- Hard discs, in 2-dimensions are 
locally jammed if they have at least 3 contacts that are 
not all in the same semicircle. By restricting the chan- 
nel radius to values oi h < •\/3/4 we ensure particles 
can only interact with the two discs on either side so a 
jammed particle has two disc-disc contacts and one disc- 
wall contact. As a result, there are only two types of 
particle arrangement that lead to locally jammed parti- 
cles (see Fig. ^) . The most dense arrangement requires 
the contacting neighbors to form a "V" with the central 
disc so that a series of densely packed discs forms a zig- 
zag structure. A defect, or locally less dense packing has 
two particles on the same side of the tube. The defect 
behaves like a vacancy in that a particle in the defect 
can hop into the empty "lattice" site but it is not ac- 
tually possible to add another particle to the system at 
the defect unless h > 1^3/4. Later, we will see that the 
presence of these defect sites plays an integral role in the 
glassy behavior of the fluid at high densities. The collec- 
tively jammed configurations of the entire svstemfTsjl. or 
inherent structures are made up of combinations of these 
local packing arrangements, excluding those that contain 
neighboring defects because, as Fig. ^ shows, the cen- 
tral disc in a divacancy is free to move vertically which 
would allow the packing to unjam. 

To count the number of packings with a density zq we 
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FIG. 1: (a) Stable packing configurations. The most dense 
arrangements are white. Particles in a defect appear grey, (b) 
The loose particle in a divacancy allows the system to unjam 
and eliminate the two defects. The "0"s and "1" denote the 
bond assignments for the configurations. 



note that, by drawing a "bond" between a particle and its 
neighbor on the right and writing a "1" corresponding to 
a defect bond (i.e. when two discs are jammed together 
on the same side of the channel) and a "0" otherwise, 
we can develop a lattice gas description of a jammed 
configuration, as shown in Fig. ^ The total number of 
bonds must equal the number of particles N. If M is the 
number of defect bonds we can divide a configuration 
into blocks of zeros and ones with Mqi 01 boundaries. 
The total number of possible configurations is given by 

^ M\{N-My. 

" [M - (Moi/2)]! [N - M - (Moi/2)]! [(Afoi/2)!]' 

(1) 

If we allow all configurations, including unstable arrange- 
ments such as "11" and "111" etc, then Eq.^ becomes 
the ising model expression Ng{zo) = N\/M\{N — A/)!. 
To count only the stable packings, we remove all those 
configurations including neighboring defects, or "11" ar- 
rangements, by enforcing Afoi — 2Af so that Eq. ^ re- 
duces to 

(iV- M)! 



M\{N -2M)\ 



(2) 



The density of a jammed state is given by zq = 

[(l-6')Vl-/i2 + 6']~\ where 9 = M/N is the mole 
fraction of defects, and the configurational entropy is 
SjNk = {l/N) \nNg{zo) = (1-0) ln(l-6')-6'ln6l- (1- 
29) ln(l — 29). The distribution of inherent structures is 
binomial with a single close packed structure with no de- 
fects and a single least packed structure with 9 = Q.b and 
maximum number of inherent structures occurring with 
9 = 1/2-^/5/10. The same binomial distribution of in- 
herent structures is found for a mixture of non-additive 
one-dimensional hard rodsQ- 

To complete the density landscape picture, we need 
the contribution to the entropy associated with the con- 
figurations that map to each inherent structure. The 
partition function of an individual basin could be calcu- 
lated numerically using constrained simulations but, 
for our immediate purposes, it is sufficient to make use 
of the simple double line model originally introduced by 
Wojciechowski et al 19] which maps the model onto a 
mixture of 1-dimensional hard rods that interact with 
their neighbors with diameter a if both discs are on the 
same side of the channel or cr(l — h^) if they are not. The 
resulting 1-dimensional approximation for the partition 
function of a basin with an inherent structure of density 
Zq is Qg{z,ZQ) = (A^L^/7V!)(zo - z)/zq where A is the 
thermal DeBrogiie wavelength of a disc. The equilibrium 
properties of the fiuid are given by maximising the to- 
tal entropy so that Sf{z)/kN — IniVg(zQ) -I- lnQg{z, Zq), 
where Zq denotes the value of the limiting density that 
satisfies {dSf /Ozq)^ — 0. 

A key feature of the model is that the number and 
position of the defects indentify which basin the system 
is visiting. One method of counting the defect concen- 
tration in the fluid is to periodically interrupt the sim- 
ulation and compress the configuration to its inherent 
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structure. However, due to the nature of the model, it 
is actually possible to determine the lattice model de- 
scription of the inherent structure from an instantaneous 
configuration by considering the relative positions of the 
neighboring particles and using the triangular constraint 
used by SpeedyQ to examine configurations of bulk disc 
systems. For example, if a particle is below the line drawn 
between the centres of its neighboring discs, the particle 
will pack against the bottom wall, otherwise it becomes 
jammed against the upper wall. Fig. ^ shows the aver- 
age concentration of defects obtained from a simulation 
of N = 5000 particles as a function of density compared 
to that predicted by the model. At low densities, the 
system visits the set of basins that maximizes its con- 
figurational entropy but, at higher densities, the defects 
are eliminated in order to create more free volume. Due 
to the binomial distribution of inherent structures, the 
ideal glass transition for the model only occurs at close 
packing when the number of defects goes to zero. Some 
di vacancies were also observed in the simulations. For a 
system with h = 0.866, at 2; = 0.11, about 3.6% of the 
defects appear in divacancies and this concentration de- 
creases below 0.5% by z = 0.5. The insert shows that 
the model improves as the channel narrows because the 
1-dimensional partition function becomes a better ap- 
proximation. 

Recent simulation studies of supercooled fiuids have 
focused on relating structural relaxation to the organisa- 
tion of the PEL into superstructures called metabasins. 
However, the notion of a metabasin is difficult to de- 
fine. The system is usually considered to have escaped 
from one metabasin to another during s simulation when 
a considerable change in the potential energy of the 
quenched structure correlates with changes in the par- 
ticle positions 20] . The relaxation of the system is then 
studied in terms of the distribution of waiting times for 
hops between metabasins |23- In the present model we 
can group all the basins belonging to inherent structures 
with the same number of defects into a single metabasin. 
The system moves between basins in the same metabasin 
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FIG. 2: measured from simulation (points) compared to the 
theory (lines) for a channel of width h = 0.866, as a function 
of z The insert: h = 0.5. 
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FIG. 3: Distribution of waiting times. From the right z — 
1.306, 1.188, 1.069, 0.950, 0.831, 0.712, 0.593, 0.475, 0.356. The 
lines represent best fits to the expression for P{t) in the text. 
The inset shows the best fit parameter /3 as a function of 
density. 



when a particle hops to the vacant site in a defect. The 
transition between metabasins occurs through configu- 
rations containing a divacancy. Fig. shows that the 
translation of the central disc results in the elimination of 
two defects. The creation of defects occurs via the same 
transition state. 

We can also directly identify the configurational ele- 
ments that lead to dynamic heterogeneity and local rules 
for dynamic facilitation of defect hoping and defect elim- 
ination events. As the density of the fiuid increases, the 
particles become caged by their neighbors and so they can 
only rattle around their local lattice site. However, the 
defects act as local regions of "excitation" and particles 
located in a defect can also hop to the vacant site. This 
is a very simple example of facilitation in which the local 
dynamics rules arise out of local packing constraints. As 
yet, we have not identified a configurational clement that 
leads to the creation of defects. 

Fig. 131 shows the distribution of waiting times for tran- 
sitions between metabasins on the PEL, i.e. when the 
number of defects in the system changes. At high den- 
sities, where the caging of particles means that hopping 
at a defect is the only way to move between basins, we 
see long relaxation times that are associated with the 
need for defects to diffuse together and the distribution 
is well described by the facilitated dynamics model with 
P{t) cx (t/trei)^ exp [—{t/trei)^], whcrc trei and f3 are 
both best fit parameters. The insert shows that the ex- 
ponent for the stretched exponential becomes linearly de- 
pendent on density at high densities and in fact, the best 
fit line suggests that f3 would go to zero at close pack- 
ing. At low densities, where there are no barriers between 
basins on the landscape and no local packing constraints 
to require hopping between basins to be facilitated, the 
expression for P{t) is no longer a good fit for the waiting 
times distribution. 

The simplicity of our model allows us to examine 
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FIG. 4: Exchange times and Persistence times {Pp{t)) 

distribution functions, (a) tpx{t)a.t z = 0.712 (b) tp^it) at 
z = 1.425. (c) tpp(t) ak. z = 0.712 (d) tpp{t) at z ^ 1.425. 
Heavy lines denote the full distribution. Thin lines and dotted 
lines denote excitation (0 l)and de-excitation (1 ^ 0) 
events respectively. 

some of the properties of the exchange and persistence 
times ^3 that help characterize heterogeneous dynam- 
ics. Fig. ^ shows the distributions of exchange times, 
defined as the time between successive events (either ex- 
citation ^ 1 or a de-excitation 1 ^ 0) at a particular 
position in the fluid, and the persistence times, which 
measure the time a local region takes before it is either 
excited or de-excited for the first time. We used systems 
of ~ 1006*"^. Breaking the distributions into their 
components shows that the slow relaxation at high den- 
sities is dominated by the presence of well packed regions 
waiting to be excited by the appearance of a defect which 
is rare and must diffuse through the system. 

The exchange time distributions exhibit a rich struc- 
ture at high densities that is associated with processes 
of defect creation, elimination and hopping that we will 
not describe in this letter. However, from the exchange 
times, we can extract an average hopping time associated 



with the time a defect takes to jump to its neighboring 
site by considering just the de-excitations but exclud- 
ing any elimination events. Fig El shows that the av- 
erage hopping time for a defect fits a free volume law, 
r = A exp[B z / {zq — z)] above z — 0.8 where B Ri 1. This 
can be expected on physical grounds because, in order 
to hop from one site to the next, the cage must open 
sufficiently to allow the particle to move between its two 
neighbors. 

In summary, we have presented a simple off-lattice 
model for which we can directly identify important con- 
figurational components and local dynamical rules, that 
lead to heterogeneous dynamics and facilitation, and are 
directly related to local packing constraints. This is 
consistent with the notion suggested by the propensity 
of particles in more complex systems'll^, that dynamic 
heterogeneity arises from a configurational component. 
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FIG. 5: Free volume plot for the average hopping times. 



We have shown that our model displays glassy dynam- 
ics purely in terms of the landscape paradigm through 
its distribution of metabasin waiting times. However, we 
have also shown that the dynamics of our model can be 
described in terms of dynamic facilitation. 
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